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Abstract 



Using recent values of the QCD (non-) perturbative parameters given in Table 1 and an estimate of the N3L0 QCD perturbative 
contributions based on the geometric growth of the PT series, we re-use QCD spectral sum rules (QSSR) known to N2L0 PT series 
" and including all dimension- six NP condensate contributions in the full QCD theory, for improving the existing estimates of rhc^b 
^-H and fD(,),B(,) from the open charm and beauty systems. We especially study the effects of the subtraction point on "different QSSR 
data" and use (for the first time) the Renormalization Group Invariant (RGI) scale independent quark masses in the analysis. The 
(N estimates [rigourous model-independent upper bounds] reported in Table 8: fo/fn = 1.56(5)[< 1.68(1)], /fi//;^ = 1.58(5)[< 1.80(3)] 
^LiSLYidfoJfK = 1.58(4)[< 1.63(1)], /b J /k = 1.50(3)[< 1.61(3.5)], which improve previous QSSR estimates, are in perfect agreement 
D "(in values and precisions) with some of the experimental data on f^D, and on recent lattice simulations within dynamical quarks. 
These remarkable agreements confirm both the success of the QSSR semi-approximate approach based on the OPE in terms of the 
quark and gluon condensates and of the Minimal Duality Ansatz (MDA) for parametrizing the hadronic spectral function which we 
have tested from the complete data of the J/i// and T systems. The values of the running quark masses mdmc) = 1286(66) MeV and 
rfihimh) = 4236(69) MeV from Md,b are in good agreement though less accurate than the ones from recent J/i// and T sum rules. 
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1. Introduction and a short historical review 

The (pseudo)scalar meson decay constants fp are of prime 
interests for understanding the realizations of chiral symmetry 
in QCD. In addition to the well-known values of fj^ = (130.4(2) 
MeV and //^ = 1 56. 1 (9) MeV [2] which control the Hght flavour 
.chiral symmetries, it is also desirable to extract the ones of 
the heavy-light charm and bottom quark systems with high- 
accuracy. These decay constants are normalized through the 
matrix element: 



: {ou^Qixm = fpMi , 

where: 



(1) 



(2) 



is the local heavy-light pseudoscalar current; q = d, s; Q = 
c,b; P = D(s), ^(^)and where fp is related to the leptonic width: 



(3) 



where m/ is the lepton mass and \VQq\ the CKM mixing angle. 
Besides some earlier attempts based on non-relativistic poten- 
tial models to extract these quantities (which are however not 
applicable for the heavy-light systems), the first bounds on fo 
and /b from QCD spectral sum rules (QSSR) [3] ^ were derived 



*Some results of this work have been presented at the 16th QCD Interna- 
tional Conference (QCD 12) , Montpellier, 2-6th July 2012 [1]. 
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by NSV2Z [8], which have been improved four years later in 
[9-1 1]. Since then, but long time before the lattice results, dif- 
ferent QSSR papers have been published in the literature for 
estimating fo^B ^- These results look, at first sight, in disagree- 
ment among each others and some of them, claimed the ob- 
servation of the scaling /p ~ 1 / ^fM^ expected in the large 
Mp limit [12]. These diff'erent papers have been scrutinized 
in [6, 13], where Narison found that the apparent discrepancies 
between the diff'erent results can be solved if one applies care- 
fully the stability criteria (also called sum rule window) of the 
results versus the external QSSR Laplace/Moments sum rules 
variables and continuum threshold tc. In this way, and for given 
values of m^b, he obtained the values: 



fn - (1.31 ± 0.12)A , fp - (1.6 ± 0.1)A . 



(4) 



which are independent of the forms of the sum rules used. 
However, the result has been quite surprising as it indicates a 
large violation of the heavy quark symmetry scaling predic- 
tions,where 1 /Mq corrections have been estimated in [14]. This 
"unexpected result" has been confirmed few years later by lat- 
tice calculations [15]. Since then, some progresses have been 
done for improving the QCD expression of the 2-point corre- 
lator. It starts from a confirmation of the SVZ original ex- 
pression of the LO perturbative and non-perturbative contribu- 
tions. Then, Broadhurst and Generalis [10, 16] have provided 
the complete PT as NLO including light quark mass correc- 
tions. It has been completed by the PT a^^ N2L0 corrections 



^ For reviews and more complete references, see e.g: [5, 6]. 
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of Chetyrkin and Steinhauser [17] in the case of one heavy and 
one massless quarks. This result has been completed by the in- 
clusion of the NP contributions up to dimension-six [11] and of 
the light quark mass corrections to LO by [1 1, 18]. All of these 
previous QCD expressions have been given in terms of the on- 
shell quark mass. In [19], Narison has used (for the first time) 
the running c, b quark masses in the QSSR analysis, by using its 
known relation with the on-shell mass known at present to NLO 
[20-22], N2L0 [10, 16] and N3L0 [24] where it has been no- 
ticed that the QSSR PT expressions converge faster. It has also 
been noticed that the values of fD,B are very sensitive to the 
value of ntcb motivating him to extract mc,b (for the first time) 
from the known values of Md and Mb. Recent analysis, includ- 
ing the corrections have been presented in the literature, in 
the full theory where the running MS mass has been used [25- 
27] and in HQET [28] where the radiative corrections are large 
due to the uses of the on-shell mass ^. 

In the following, we shall present analysis based on the full 
QCD theory where we use as inputs the most recent values of 
the (non-)perturbative QCD parameters given in Table 1 . We 
assume the geometric growth of the PT series [30] as a dual 
to the eff'ect of a l/q^ term [31, 32] for an estimate of the 
N3L0 perturbative contributions. We shall also study systemat- 
ically the eff'ect of the substraction points on each "QSSR data" 
and use (for the first time) in the analysis, the Renormalization 
Group Invariant (RGI) s, c, b quark masses introduced by [33] 
and which are scale and (massless) scheme independent. 

2. QCD spectral sum rules (QSSR) 

• The Laplace sum rules (LSR) 

We shall be concerned with the two-point correlator : 



J d'x e''^-'{0\TJlQ{x)JlQ{0)^\0) , (5) 

where JqQ(x) is the local current defined in Eq. (2). The as- 
sociated Laplace sum rules (LSR) JLqQij) and its ratio "RqQir) 
read [3]^: 



f 

J(m, 



,2 dt t e~^^-\mil/^-Jt, u) 

J{mq+MQ)^ n ^qQ^ 



(6) 



(7) 



where fi is the subtraction point which appears in the approxi- 
mate QCD series when radiative corrections are included. The 
ratio of sum rules *R^g(T,yu) is useful, as it is equal to the res- 
onance mass squared, in the Minimal Duality Ansatz (MDA) 
parametrization of the spectral function: 

- Imj/r^^(0 ^ fpM^pS(t - Ml) -h "QCD cont."^(^ - tc)X^) 



^We plan to analyze the HQET sum rules [14, 29] in a separate publication. 

"^Radiative corrections to the exponential sum rules have been first derived 
in [34], where it has been noticed that the PT series has the property of an 
Inverse Laplace transform. 



where fp is the decay constant defined in Eq. (1) and the higher 
states contributions are smeared by the "QCD continuum" com- 
ing from the discontinuity of the QCD diagrams and starting 
from a constant threshold tc. 

• The = moment sum rules (MSR) 

We shall also use for the 5-meson, the moments obtained after 
deriving n -h 1 -times the two-point function and evaluated at 
= [3], where an expansion in terms of the on-shell mass 
Mh can be used. They read: 



and the associated ratio: 

s. 



(10) 
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Figure 1 : a) r-behaviour of the ratio of Xjf /Xj"' for = M^(25)-0. 15 GeV. The red 
dashed curve corresponds to the strict equality for all values of r. ; b) the same as a) but for 



• Test of the Minimal Duality Ansatz (MDA) from J I if/ and T 

We have checked explicitly in [6] that the MDA presented in 
Eq. (8), when applied to the p-meson reproduces within 15% 
accuracy the ratio measured from the total cross-section 
e^e~ 1=1 hadrons data (Fig. 5.6 of [6]). In the case of 
charmonium, we have also compared from 1^^^ with the 
one from complete data and find a remarkable agreement for 
higher n> A values (Fig. 9.1 of [6]), indicating that for heavy 
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Figure 2: The same as in Fig. 1 but for the Z?-quark and for ^[^c = My(25)-0.15 GeV. 



quark systems the role of the QCD continuum will be smaller 
than in the case of light quarks and the exponential weight or 
high number of derivatives suppresses efficiently the QCD con- 
tinuum contribution but enhances the one of the lowest ground 
state in the spectral integral. We redo the test done for charmo- 
nium in Fig. 9.1 of [6] and analyze the bottomium channel for 
the LSR and MSR. We show in Fig. (la) the r-behaviour of the 
ratio of Xf-^^ normalized to X^"^^ where we have used the sim- 
plest QCD continuum expression for massless quarks to order 
al from the threshold tr ^ : 



QCD cont. = l-\-as-\- l.Sas^ - 12.07 as\ 



(11) 



We show in Fig. (lb) the r-behaviour of M^, where the contin- 
uous (oliva) curve corresponds to ^/t^ ^ ^(^(25) - 0.15 GeV. We 
show a similar analysis for the bottomium sum rules in Fig. (2) 
for the LSR and in Fig. (3) for the MSR where we have taken 
^/t^ ~ Mr (28) - 0.15 GeV. One can see that the MDA, with a 
value of ^/t^ around the value of the 1st radial excitation mass, 
describes quite well the complete data in the region of r and n 
where the corresponding sum rules present r or ^ stability [35]: 



- (1.3 
^ (0.2 



^ 1.4) GeV-^ 
^ 0.4) GeV-^ 



^ (5 



(12) 



as we shall see later on. This good description of the data by 
the MDA shows the efficient role of the exponential weight 
or high number of derivatives for suppressing the higher mass 



^We have checked that the spectral function including complete mass cor- 
rections give the same results. 



States and QCD continuum contribution in the analysis. This 
nice feature prevents the introduction of some more involved 
models bringing new parameters in the analysis where some 
of them cannot be understood from QCD 1st principles. More- 
over, MDA has been also used in [36] (called Minimal Hadronic 
Ansatz in this paper) in the context of large Nc QCD, where 
the restriction of an infinite set of large Nc narrow states to a 
Minimal Hadronic Ansatz which is needed to satisfy the lead- 
ing short- and long-distance behaviours o the relevant Green's 
functions, provides a very good approximation to the observ- 
ables one compute. 
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Figure 3: The same as in Fig. 2 but for the = moment of the Z^-quark versus the 
number of derivatives n. 
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Figure 4: a) r-behaviour of "^(r) normalized to the ground state energy £"0 for the har- 
monic oscillator 2 and 4 indicate the number of terms in the approximate series. 
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• Optimal results from stability criteria 

Using the theoretical expressions of or M^^^^^, and 

parametrizing its experimental side or M^j^^^^ by the 
MDA in Eq. (8), one can extract the decay constant fp and 
the RGI quark mass rriQ. In principle the equality Xj^ = 
should be satisfied for any values of the external (unphysical) 
set of variables (r, tc), if one knows exactly Xj^ and Xj^. Un- 
like the harmonic oscillator, this is not the case. Using the ratio 
of moments "Rjq for the harmonic oscillator as a function of the 
imaginary time variable r, where one knows the exact and ap- 
proximate results, one can find [37] that the exact energy £"0 of 
the ground state can be approached from above by the approx- 
imate series (see Fig. 4). At the minimum or inflexion point 
(stability) of the curves, one has a ground state dominance. For 
small time (large Q^), all level contributes, while for large time 
(small Q^) the series breakdown. We shall apply this stability 
criterion inspired from quantum mechanics in our analysis. 
In principle, the continuum threshold ^/t^ in Eq. (8) is a free 
parameter, though one expects its value to be around the mass 
of the 1st radial excitation because the QCD spectral function 
is supposed to smear all the higher state contributions in the 
spectral integral as explicitly shown previously in Section 2. 
In order to avoid the model-dependence on the results, Refs. 
[5, 6, 13, 14, 19, 25] have considered the conservative range of 
-values where one starts to have r- or ^-stability until which 
one reaches a stability where the contribution of the lowest 
ground state to the spectral integral completely dominates. For 
the D and B mesons, this range is [5, 6, 13, 14, 19, 25]: 



t^ ^ (5.5 ~ 9.5) GeV^ ^ (33 ~ 45) GeV^ 



3. The QCD input parameters 



(13) 



The QCD parameters which shall appear in the following 
analysis will be the strange, charm and bottom quark masses 
^s,c,b (we shall neglect the light quark masses q = u, J), the 
light quark condensate (qq), the gluon condensates {g^G^} = 
{g'Gl.G^J) and {g^G^) ^ {g' fatcG^^^G'^^G^^X the mixed con- 
densate (qgcrGq) = {qgcr^''(Aa/2)G^^q} = M^iqq) and the four- 
quark condensate p{qq)^, where p ^ 2 indicates the deviation 
from the four-quark vacuum saturation. Their values are given 
in Table 1 and we shall work with the running light quark pa- 
rameters known to order a] [5, 6, 38]. They read: 

rhq^Qir) = rhq^Q {-pxas)~^'^' X C{as) 
{qq){T) = -il\{-p,asflP'IC{as) 
(qgaGqKr) = -Mo'//^ (-^ia,)i/^^VC(a,) , (14) 

where fii = -(1/2)(11 - 2n//3) is the first coefficient of the /3 
function for rif flavours; as = as(T)/7T; ihq^Q is the RGI quark 
mass, jlq is spontaneous RGI light quark condensate [33]. The 
QCD correction factor C(a^) in the previous expressions is nu- 
merically: 

C{as) = 1 + 0.8951a, 1.3715aJ ... : n/ = 3 , 

= 1 1.1755a, 1.5008aJ ... : rif = 5 , (15) 



which shows a good convergence. We shall use: 
as(Mr) = 0.325(8) =^ as(Mz) = 0.1192(10) 



(16) 



from T-decays [39, 40], which agree perfectly with the world 
average 20 12 [41, 42]: 



a,(Mz) = 0.1184(7). 



(17) 



We shall also use the value of the running strange quark mass 
obtained in [43] ^ given in Table 1. The value of the running 
(qq) condensate is deduced from the value of (m^ + mj)(2) = 
(7.9+0.6) MeV obtained in [43] and the wefl-known GMOR re- 
lation: (mu + md){uu + dd) = -m^f^. The values of the running 
MS mass mQ(MQ) recently obtained in Ref. [35] from char- 
monium and bottomium sum rules, will also be used^. Their 
average is given in Table 1. From which, we deduce the RGI 
invariant heavy quark masses to order a^, in units of MeV: 



rhc = 1467(14) , 



m = 7292(14) . 



(18) 



For the light quarks, we shall use the value of the RGI mass and 
spontaneous mass to order a, for consistency with the known a, 
and (qq) condensate corrections of the two-point correlator. 
They read, in units of MeV: 



ms = 128(7) , 



flq= 251(6). 



(19) 



Table 1: QCD input parameters. 



Parameters 



Values 



Ref. 



a, (Mr) 

m,(2) 

riicimc) 

rhbimb) _ 

\{uu + ddyi\2) 

{ss)l{dd) 

Ml 

{asG^) 
ig'G') 
p{qq)^ 



0.325(8) 
96.1(4.8) MeV 
1261(12) MeV 
4177(11) MeV 
-(275.7 + 6.6) MeV 
0.74(3) 

(0.8 + 0.2) GeV^ 
(7 + 1) X 10-2 GeV^ 
(8.2 + 1.0) GeV^ X (a.G^) 
(4.5 + 0.3) X 10-4 GeV^ 



[39-41] 
average [43] 
average [35] 
average [35] 
[5,43] 
[5, 43, 47] 
[48-50] 

[35,37, 39,51-56] 
[35] 

[39, 48,51] 



4. QCD expressions of the sum rules 
• The LSR 

To order a^, the QCD theoretical side of the sum rule reads, in 
terms of the on-shell heavy quark mass Mq and for md = 0: 



jl+2a,[l + (l-z)|lnv2T+^j 



^This value agrees and improves previous sum rules results [44]. 
^These values agree and improve previous sum rules results [3-6, 45, 46]. 



4 



-2a,r(0,z)M^^j maidd). 



-TC-4^il-^)MQMl{dd} 



+ 2 



z^ydjd) 



5L(l2-3z-z^)z-9 + Uz+^z 



where: 



PT 



1 



(f) ] 

21607r2J ' 



3t(l-xy\l + -asf(x)] + atR2s 



(20) 



(21) 



with: z = MqT; x = M^qIv, Us = asl^\ L = InipiMgr) + : 
Ye = 0.577215...; yu is an arbitrary subtraction point; 7^2^ is the 
Qf^-term obtained semi -analytically in [17] and is available as a 
Mathematica package program Rvs.m. Neglecting mj, the PT 
NLO terms read [10]: 

9 

fix) = - + 2Li2(x) + log X log( 1 - x) 
-^log(l/x-l)-log(l-x) + 



+ ln|^j (a, + (8.8472- 0.361 ln/)a^) 

+ ln2|-^j (1.7917 -0.0833fz/) a? + ... , (24) 



for ni light flavours, one can express ah terms of the previous 
sum rules with the running mass mqiji). It is clear that, for 
some non-perturbative terms which are known to leading order 
of perturbation theory, one can use either the running or the 
pole mass. However, we shall see that this distinction does not 
aff'ect, in a visible way, the present result, within the accuracy of 
our estimate, as the non-perturbative contributions are relatively 
small though vital in the analysis. 

• The MSR 

The moments read for rud = 0: 



•(n) 
db 



1 



+{n + l){n + 2)——M^{dd) 



xlog(l/x - 1) - (x/(l - x))\ogx . 



(22) 



4M^ 

-{n + l){n + 2){n + 9) 



-{n + 3){5rf + 9n + 1) 



1 {djd) 
Ml 36 

ml 28807^2 



The contribution up to the J = 4 gluon condensate and up to J = 
6 quark condensates have been obtained originally by NSV2Z 
[8]. The contribution of the J = 6 (g^fabcG^) and (f) gluon 
condensates have been deduced from the expressions given by 
[11] (Eqs. II.4.28 and Table II.8) where: 

(djd) ^ 0gy^D^G,y^d) = gHdy,^dJ]qy,^q) 



16 - 

^ -—{7ias)p{dd) , 

- -^{nasfp{dd)\ 



\2 



V q 



(23) 



after the use of the equation of motion, p ^ (2 + 0.2) measures 
the deviation from the vacuum saturation estimate of the J = 6 
quark condensates [39, 48, 51]. 

The as correction to (Jj), in the M^" -scheme, comes from [26], 
where the running heavy quark mass mq enters into this expres- 
sion. Using the known relation between the running friQifi) and 
on-shell mass Mq in the M5 -scheme to order or^ [20-24]: 



Mq = mQijj) 



1 + -a, + (16.2163 - 1.0414n/)aJ 



where: 
S 



-(20n^ + 186n^ + 337n + 117) + 



-5(n + 2) 



^4(^^ + 7^+ 12) + 
3^3(^ + 3)- 12^2- 



(n^+ 10n + 9)ln| — 
A^ 



1 (/) ] 
Ml 216071^ j 



, (25) 



n ^ 



(26) 



5. Estimates of fp and mQSLtfi = r from LSR 

After inspection, one finds that fp and the RGI mass niQ 
can only be simultaneously determined from XjQ(T,yu) and 
!Rjg(T,yu) evaluated at /i = t"^/^. For other values of /i, only 
£,j^(t,ij) present r stability at reasonable values of r < 1.2 
GeV~2, which is not 'Rjci'^, ji). This particular value of 11 = I ^ 
is also interesting because the subtraction scale moves with the 
sum rule variable r in the analysis. 
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Figure 5: a) r-behaviour of fo from Xjc for different values of tc, for a given value of the 
subtraction point fi - t~^^^ and for the - 1467 MeV as given in Eq. (18); b) the same as a) 
but for Md from "/^j^. 



• Analysis of the r- stabilities ofJlj^ and *Rj^ 

Using the central values of the QCD input parameters in Table 
1 and in Eqs. (16), (18) and (19), one can show in Fig. 5 the 
influences of r and tc on the value of fo and Md for a given 



value of the subtraction point // 
for //) is reached for: 



-1/2 



where, the r-stability 



(0.8 ~ 1.2) GeV-2 , t^ ^ (5.3 ~ 6.5) GeV^ (27) 



When extracting the RGI mass rhc from *Rj^ by requiring that it 
reproduces the experimental mass squared M^, one can notice 
in Fig. (5) that, unlike fo, Md present r-stabilty for larger range 
of ^c-values^: 



^ (5.3 ~ 9.5) GeV^ 



(28) 



• Analysis of the convergence of the QCD series 

We study the convergence of the QCD series in the case of 
the charm quark at a such low value of the subtraction point 
[1 = T"^/^ and taking tc = 6 GeV^ . We work in the M^" -scheme 
as we know from previous works [19] that the PT series con- 
verge better than using the on- shell subtraction. In so doing, we 
estimate the a] N3L0 contribution using a geometric PT se- 
ries as advocated in [30] which is dual to the efl'ect of the 1 /q^ 
term when large order PT series are resummed. We show the 
T-behaviour of fo in Fig. (6). One can notice that, all correc- 
tions act in a positive way. The prediction increases by about 
17% from LO to NLO and another 14% from NLO to N2L0 
but remains unafl'ected by the inclusion of the N3L0 contribu- 
tion estimated above. These features indicate that the PT series 



^The existence of r-stability at values of tc below 5.3 GeV^ depends on 
the heavy quark mass value and disappears when we require the sum rule to 
reproduce the value of Md, such that we shall not consider a such region. 



converge quite well at this low scale, while the size of each PT 
corrections are reasonably small and will be even smaller for 
higher values of the subtraction point ji and for the 5-meson. 
Therefore, a confirmation of this N3L0 estimate requires an 
explicit evaluation of this contribution. 

As far as the non-perturbative contributions are concerned, their 
eff'ects are relatively small. 
0.24 1 
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Figure 6: r-behaviour of fo from LSR for GeV^, for rhc = 1467 MeV, for a given 
value of the subtraction point fi = t~^^^ GeV and for different truncations of the QCD PT 
series, where the estimated N3L0 contribution is small indicating a good convergence of 
the series; b) the same as a) but for Mg. 



• QCD and systematic error estimates 

Using the previous QCD input parameters and their correspond- 
ing errors, we deduce the diff'erent errors on fp and niQ given 
in Table 2, where the optimal results have been taken at the r- 
stability region: 



(0.8 ~ 1.2)GeV- 



(29) 



In order to take into account the systematics of the approach, 
which is gouvernedby the choice of t^ we take the central value 
of fo in Table 2 as coming from an arithmetic average of its val- 
ues from the diff'erent tc given in the legend of Fig. (5) inside 
the range given by Eq. (27). The existence of r-stability at 
values of tc below 5.3 GeV^ depends on the heavy quark mass 
value and disappears when we require the sum rule "Rir) to re- 
produce the value of Md, such that we shall not consider a such 
region. In this way, we obtain the error induced by the range 
of ^c-values given in Table 2 which we consider as systematic 
errors of the approach. As mentioned earlier, we do not try to 
introduce any model-dependent new parameters for improving 
the MDA for the continuum contribution because we want to 
have a stronger constraint on the results. In fact, in addition 
to the test performed in Section 2, we also want to test (a pos- 
teriori) the degree of accuracy of this parametrization for the 
heavy-light systems by confronting the results obtained in this 
paper with the some known data on fp or from lattice simula- 
tions. Indeed, diff'erent tests of this MDA model are important 
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as it is widely used in the literature for predicting not yet mea- 
sured masses of new hadrons like four-quark, molecules [57] 
and hybrid [58] states. However, we do not also try to fix more 
precisely tc by e.g. using Finite Energy Sum Rule [53] like did 
Refs. [59] as we want to have a more conservative result. 

• Results for fo and rhc 

Considering the common range of -values for fo and Md 
given in Eq. (27), we obtain the results quoted in Table 2 which 
come from an arithmetic average of optimal values obtained at 
diff'erent tc values in Eq. (27) ^: 

fo = 204(11) MeV, 

rhc = 1490(77) => m,(m,) = 1286(66) MeV . (30) 

which we consider as improvement of the result obtained from 
the same sum rule and at the same subtraction point by [19] 

fo = 205(20) MeV , mdmc) = 1 100(40) MeV . (31) 

The smaller errors in the present analysis, come from more pre- 
cise input parameters, more complete NP-corrections included 
into the OPE and more constrained range of ^^-values.The value 
obtained in Eq. (30) also agrees within errors with the ac- 
curate determination from charmonium systems quoted in Ta- 
ble 1 though less accurate. The main sources of errors from the 
present determination can be found in Table 2. One can notice 
that the contributions of the J = 6 condensates are negligible 
for fo (less than 0.3 MeV) and small for ({ddf and <G^) 
contribute respectively to 17 and 6 MeV). 

Table 2: Central values and corresponding errors for fp and mq in units of MeV 
from the LSR at the subtraction point /u = t~^^^. We have used mg in Eq. (18) 
for getting fp. The +(resp. -) sign means that the values of fp, mq increase 
(resp. decrease) when the input increases (resp. decreases). The relative change 
of sign from c to bin some errors is due to the effects of rm^ appearing the OPE. 
Notice that the error in (G^> also affects the (G^> contribution. The Total error 
comes from a quadratic sum. 





Value 


tc 


a, 




niQ (dd) (G2> 


Ml 


(dd)^ (G^) Total 


fo 


204 


+4 


-9 


+3 


-2 +3.5 +0.5 


-0.5 


-0.01 +0.03 11 


fB 


201 


+7 


-10 


+ 1 


-2 +1.9 +0.05 


-0.25 


-0.00 +0.00 13 


rhc 


1457 


-44 


-64 


-24 


+22 +5 


-38 


+ 1.5 -0.8 93 


rhb 


7272 


-150 


-114 


-14 


+20 +5 


-39 


-13 -14 195 



• Extension of the analysis to fB and rhb 

We extend the previous analysis to the case of the /?-quark. We 
show in Fig. (7) the r-behaviours of /g and Mb for different 
values of tc. One can see, that in this channel, r- stability for /g 



^ Using the larger range of -values, we would have obtained a slightly dif- 
ferent value: rhc - 1492(102)^^(82)^cJ MeV, where the errors come respec- 
tively from the choice of tc and QCD parameters given in details in Table 2. 

^*^An extended discussion about the value of fo at different subtraction points 
will be done in the next section. 
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Figure 7: a) r-behaviour of from Xj^ for different values of tc, for a given value of the 
subtraction point = r"^/^ and for % = 7292 MeV as given in Eq. (18); b) the same as a) 
but for Mb from "^j^. 



is reached at ^ . 

tI ^ (0.2 ~ 0.26) GeV"^ , ^ (33 ~ 35) GeV^ (32) 

while, like in the case of Md, Mb stabilizes for a larger range 
of values 

^ (33 ~ 45) GeV^ . (33) 

We show in Fig. (8) the predicted values of /g and Mb for a 
given value of fhij and for different truncations of the PT QCD 
series. Using the same procedure as in the charm quark case 
and considering the range of tc in Eq. (32), where the central 
values of /g and rhh, in units of MeV in Table 2 comes from an 
arithmetic average of different optimal values in the range of tc 
in Eq. (32), we deduce the estimate in units of MeV: 

fB = 201(13) 

m = 7272(195) =^ m/,(m/,) = 4164(112), (34) 

which we again consider as improvement of the result from 
[19]: 

fB = 203(23) MeV, m^(m/,) = 4050(60) MeV , (35) 
obtained from the same sum rule. 



^^The apparent minima at r < 0.1 GeV^ obtained for lower values of tc cor- 
responds to the region where the continuum contribution to the spectral integral 
is dominant and should not be considered. 

^^Like in the case of the charm quark, we shall not consider values of tc < 
32.5 GeV^ where the r-stabilty disappears, when one requires the sum rule to 
reproduce Mb. 
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Figure 8: r-behaviour of fe from LSR for t^=33 GeV^, for % = 7292 MeV, for a given 
value of the subtraction point /i = t~^I^ GeV and for different truncations of the QCD PT 
series; b) the same as a) but for Mg. 




Figure 9: a) r-behaviour of fo from LSR for different values of tc, for a given value of 
the subtraction point jj. - \A GeV and for mc=1467 MeV; b) the same as a) but for /b, 
using = 3 GeV and %=7292 MeV 



6. Effects of the subtraction point on fj),B from LSR 

The choice of subtraction points is also one large source of er- 
rors and discrepancies in the existing literature. In order to cure 
these weak points, we extract the values of fo^B and the corre- 
sponding errors at a given value of the subtraction point ji. We 
show in Fig. (9) the r-behaviour of fo^B for given values of jd 
and ihc^b- We show in Tables 3 and 4, the results of the analysis 
including the different sources of the errors, where the typical 
sizes normalized to the values of fD,B are: 
-fo'il- 8)% from tc, (0.7 ~ 2)% from the PT contributions, 
0.5% from rUc and (0.9 ~ 1.6)% from the NP-contributions. 
- fB are: (2 ~ 4)% from tc, 4% from the PT contributions, 4% 
from mij and (0.8 ~ 1.5)% from the NP-contributions. 
We show in Figs. (10) and (14), the set of "QSSR data points " 
obtained in this way for different values of yu. 



23Q| : 

220: T I 1 T 
_ 210^ I ^ 

I 200: \ t : 

S - 1 \ - 

^ 190: I \ j \ : 

180: \ 

MO i - 

1.0 1.5 2.0 2.5 3.0 

// [GeV] 

Figure 10: Values of fo from LSR at different values of the subtraction point // and for 
mc = 1467 MeV. The filled (grey) region is the average with the corresponding averaged 
errors. The dashed horizontal lines are the values if one takes the errors from the best 
determination. 



Table 3: Central values and corresponding errors for fn in units of MeV from 
the LSR at different values of the subtraction point // in units of MeV and for 
mc=1467 MeV The +(resp. -) sign means that the values of fo increase (resp. 
decrease) when the input increases (resp. decreases). The Total error comes 
from a quadratic sum. 



1^ 


fD 


tc 


as 




Jflc 


(dd) 




Ml 


(ddf 


(G^) Total 


1.4 


204 


+ 14 


-1.3 


+4 


-1 


+3 


+ 1 


+0.6 


+0.6 


+0.6 15.0 


1.8 


204 


+ 16 


-1.2 


+2.7 


-0.9 


+2.3 


+0.3 


+0.4 


0.0 


0.0 16.5 


2.2 


203 


+ 16 


-1.0 


+2.2 


-0.7 


+2.1 


+0.3 


+0.3 


0.0 


0.0 16.3 


2.6 


203 


+ 16 


-1.1 


+ 1.5 


-1.1 


+ 1.6 


+0.3 


0.0 


-0.6 


-0.5 16.2 


3.0 


201 


+ 17 


-0.8 


+ 1.2 


-0.8 


+ 1.6 


+0.3 


+0.1 


-0.5 


-0.5 17.2 



• Final results for fo and fBfrom LSR 

Using the fact that the "physical observable" is independent of 
we average (fit horizontally) the different data points of fo 
from LSR in Tables 2 and 3 and Fig. (10) (red triangle). The 



average is represented by the horizontal band in Fig. (10). The 
narrower (grey) domain corresponds to the resulting averaged 
error, while the larger one corresponds to the case where the 
error from the most precise determination has been taken. A 
similar analysis for /g from LSR has been done using the data 
in Tables 2 and 4 and Fig. (14). We deduce from this analysis, 
the final results: 

fo = 204(6) MeV = 1.56(5)/, 
fBksR = 207(8) MeV ^ 1.59(6)/,, (36) 

where the quoted errors are the averaged errors. The previous 
errors are multiplied by about 1.8 for fo and 1.65 for /g if one 
keeps the errors from the most precise determinations. 

• Final value ofrhbfrom LSR 

In addition to the sum rule subtracted at = r"^/^, we also 
notice that the sum rule subtracted at yu = M^, where the 
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Table 4: Central values and corresponding errors for in units of MeV from 
the LSR and MSR at different values of the subtraction point ji in units of GeV 
for rhb = 7292 MeV. The +(resp. -) sign means that the values of fs increase 
(resp. decrease) when the input increases (resp. decreases). The Total error 
comes from a quadratic sum. 

/b tc as a] rrib (dd) (G^) Mg (dd)^ (G^) Total 

LSR 
3 
4 
5 
6 
7 

MSR 



196 +22 
210 +23 
213 +24 

217 +24 

218 +21 



-8.0 -1.6 -1.1 +1.9 +0.1 +0.4 0.0 0.0 23.6 

-7.6 -0.3 -1.2 +1.7 +0.1 -0.2 0.0 0.0 24.3 

-7.6 +0.4 -1.2 +1.5 +0.1 0.0 0.0 0.0 25.3 

-7.3 +0.1 -1.2 +1.6 +0.1 0.0 0.0 0.0 25.2 

-7.1 +0.5 -1.0 +1.5 +0.1 0.0 0.0 0.0 22.3 



3 


183 


+7 


-16 





-2.5 


+2 





-5 








18.4 


4 


199 


+ 10 


-22 





-3 


+3 





-9 








26.1 


5 


216 


+ 11 


-19 


+ 1 


-3 


+4 





-13 








26.0 


6 


227 


+ 17 


-21 





-4 


+3 





-17 








32.3 


7 


235 


+20 


-21 


+0.5 


-3 


+4 





-20 








35.6 



6.4 

62- 
6.0- 
5.8 ; 



5.4: 
5.2 




,00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 
T [GeV-^] 

Figure 1 1 : r behaviour of Mb from LSR for different values of tc, for m^=7292 MeV and 
at the subtraction point fi - Mb. 



Table 5: Central values and corresponding errors for rhb in units of MeV from 
LSR and MSR at different values of the subtraction point ju in units of GeV. The 
+(resp. -) sign means that the values of mt increase (resp. decrease) when the 
input increases (resp. decreases). The Total error comes from a quadratic sum. 





rhb 


tc 


Of, 




(dd) 




K 


(dd)^ 




Total 


LSR 






















Mt 


7586 


-419 


-95 


-4 


+7 


+2 


-26 








431 


MSR 






















3 


7188 


-295 


-110 


-6 


+6 


+ 1 


-174 


+5 


-0.5 


360 


4 


7360 


-301 


-102 


-6 


+5 


+ 1 


-178 


+4 


-1 


365 


5 


7490 


-306 


-99 


-4 


+8 


+ 1 


-179 


+5 





368 


6 


7598 


-310 


-99 


-4 


+9 


+ 1 


-179 


+5 





372 


7 


7686 


-312 


-97 


-4 


+9 


+ 1 


-180 


+4 


-1 


374 



a) 



0.24 
0.22 
0.20 
0.18 
0.16 
0.14 
0.12 
0.10 




t,=32 GeV^ 



8 10 12 




log (ji/Miy)-tQrm disappears in the QCD expression, presents r- 
stability [see Fig. (11)] and can then provide another estimate 
of rht. The result is given in Table 5. Taking the average of this 
result with the previous one in Table 2, we deduce in units of 
MeV: 

= 7326(178) =^ m^(m^)b^ = 4195(102), (37) 

7. 0^ = moment sum rules (MSR) for the B meson 

• Convergence of the PT series 

We show in Fig. (12) the ^-behaviours of /g and Mb for dif- 
ferent values of where one can realize a good convergence 
when the N3L0 term is included.The convergence of the PT 
series is comparable with the one of LSR shown in Fig. (8). 

• Optimal values of and rhhfrom MSR 

Using a similar procedure as for the LSR, we study, in the case 
of MSR, the ^-and stabilities of /g and rht for different values 
of the subtraction point yu. The analysis is illustrated in Fig. 
(13). The results are shown in Tables 4 and 5. One can notice 
that the sum rule does not stabilize for ji < 2 GeV, while for 
other values of yu, the range of values tc at which the ^-stability 
is reached depends on the value of the subtraction point jd and 



Figure 12: a) ^-behaviour of fs from MSR for = 32 GeV^, for rhb = 7292 MeV, for a 
given value of the subtraction point ji-A GeV and for different truncations of the QCD PT 
series; b) the same as a) but for Mb. 



are inside the range 32-42 GeV^. We show the results in Table 
4 an in Fig. (14) from which we deduce the result from the 
moments in units of MeV: 

fsWsR = 203(11), 

m\MSR = 7460(164) =^ mmb)\MSR = 4272(94) , 

(38) 



8. Final values of //>, fs and rhc^t 

As a final result of the present analysis, we take the average 
of the results from LSR for f^ and m^- This result is given in 
Eq. (30. The final results for fs and % come from the average 
of the ones from LSR and MSR shown in Figs. (14) and (15) , 
which are: 

fB = 206(7) MeV ^ 1.58(5)/,, 

rhb = 7398(121) =^ m/,(m/,) = 4236(69) MeV , (39) 
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Figure 13: a)n behaviour of /b from MSR for different values of tc, for %=7292 MeV 
and at tlie subtraction point fi = 4 GeV; b) the same as a) but for Mb . 
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Figure 14: Values of /b from LSR (red triangle) and from MSR (blue open circle) at 
different values of the subtraction point ji and for %=7292 MeV. The filled (grey) region 
is the average with the corresponding averaged errors. The dashed horizontal lines are the 
values if one takes the errors from the best determination. 




Figure 15: Values of % from LSR (red triangle) and from MSR (blue open circle) at 
different values of the subtraction point fi. Same caption as in Fig. 14. 
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Figure 16: r-behaviour of fo^ from Xjc for different values of tc, for a given value of the 
subtraction point yu = 1.4 GeV and for the - 1467 MeV as given in Eq. (18). 



expressions for N2L0 and N3L0 have been used. The non- 
perturbative contributions come from the expressions given by 
[11, 18, 26] where we have taken into account corrections of 
0(ml) for the J = 4 condensates contributions while we have 
neglected the corrections for the J = 6 condensates. We 
show in Fig. (16) the r-behaviour of fo^ for different values of 
tc at given jj = lA GeV. The results for different values of ji are 
given in Table 6 and Fig (17). 



where we have used the more precise value of m^(m^) given in 
Table 1 for getting /b. One can notice that fn - /b, confirm- 
ing previous results quoted in Eq. (4). This (almost) equality 
instead of the 1 / behaviour expected from HQET has been 
qualitatively interpreted in [60] using semi-local duality, while 
in [14] large mass corrections to the HQET lowest order ex- 
pression have been found. These results are also confirmed by 
recent lattice calculations (see Table 8). 

9. SU(3) breaking and estimates of fo, and fs, 

We extend the previous analysis for extracting fo^B, by includ- 
ing the -corrections and by talcing into account the SU(3) 
breaking of the quark condensate {ss)/{dd) given in Table 1. 

• fn, from LSR 

In so doing, we use the complete PT expression in of the 
QCD spectral function given to order ck^ by [10]. The massless 



• /b, from LSR and MSR 

In this case, we only use the PT expression to order as of the 
QCD spectral function expanded up to order which is given 
by [26] . The non-perturbative contributions are the same as in 



Table 6: Central values and corresponding errors for fo, from the LSR at dif- 
ferent values of the subtraction point // and for mc=1467 MeV. The -i-(resp. -) 
sign means that the values of fp, fhc increase (resp. decrease) when the input 
increases (resp. decreases). The Total error comes from a quadratic sum. 



)W fos c^l jfic (dd) (G^) Mq Ms (ss) Total 

264 +8.2 +2.8 +3.2 +0.2 +0.8 +0.2 +0.3 +1.2 +0.5 9.4 

1.4 247 +15 +1.0 +4.6 +0.4 +1.5 +0.7 +0.9 +1.7 +1.2 16 

1.8 236 +15 +1.1 +3.3 +1.1 +1.0 +0.0 +0.3 +1.1 +0.5 15.5 

2.2 232 +16.5 +1.2 +3.5 +1.1 +2.0 +1.4 +1.5 +2.2 +1.7 17.5 

2.6 229 +17.6 +0.3 +1.8 +0.1 +1.2 +0.2 +0.5 +1.1 +0.7 18 

3.0 226 +18.4 +1.2 +1.6 +0.2 +0.7 +0.1 +0.2 +1.0 +0.4 18.6 
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Figure 17: Values of fo^ from LSR at different values of the subtraction point // and 
for mc=l461 MeV. The filled (dark blue) region is the average with the corresponding 
averaged errors. The horizontal lines are the values if one takes the errors from the best 
determination. 
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Figure 18: Values of /b, from LSR (red triangle) and from MSR (blue open circle) at 
difi'erent values of the subtraction point fi and for %=7292 MeV. Same caption as in Fig. 
17. 



the case of fo^. We show in Fig. (19) the r-behaviour and n- 
behaviour of /g^ from LSR and MSR for different values of tc at 
given yu = 4 GeV. The results for different values of ji are given 
in Table 7 and Fig. (18). 

• Results for fo^ and /b^ 

From the previous analysis , we deduce: 



/z)^ = 246(6) MeV ^1.59(5)/^ 
/^^ = 234(5)MeV^ 1.51(4)/^ 



(40) 



which, with the help of the results in Eqs. (36) and (39) lead to: 

(41) 



^ = 1.21(4), ^ = 1.14(3). 

JD JB 



These results agree within the errors with the ones obtained by 
using the semi-analytic expressions of the correlator to order 
as [61]: 



1.15(4). 



fBs 
fB 



1.16(5), 



(42) 



with data when available [2, 62] and with recent lattice simula- 
tions (see Table 8). 
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Figure 19: a) r behaviour of /b, from LSR for different values of tc, for %=7292 MeV 
and at the subtraction point fi - 4 GeV; b) the same as a) but n behaviour of /b^ from 
MSR. 



10. Rigourous model-independent upper bounds on fD(s),B(, 



Upper bounds on fo has been originally derived by NSV2Z 
[8] and improved four years later in [9-11] and more recently 
in [27] . In this paper, we shall use LSR and the positivity of the 
continuum contributions to the spectral integral for obtaining 
the upper bounds on the decay constants. The procedure will 
be similar to the estimate done in previous sections where the 
optimal bound will be obtained at the minimum or inflexion 
point of the sum rules. In the D and Z)^ -meson channels, the 
LSR present a minimum which is well localized, while in the 
B and Bs channels, the LSR present instead an inflexion point 
which induces a new error for its localization, in addition to the 
errors induced by the QCD parameters which are the same as 
in the estimate of fD,B done in the previous sections. We show 
the results of the analysis for difl'erent values of the subtraction 
points in Fig. (20) from which we deduce the final results: 



fo < 218.4(1.4) MeV = 1.68(1)/;, 
fB < 235.3(3.8) MeV = 1.80(3)/;, 



(43) 



and: 



fo^ < 253.7(1.5)MeV^ 1.61(1)/^ 
fB, < 251.3(5.5)MeV^ 1.61(4)/^ 



(44) 



These bounds are stronger than earlier results in [9-11], while 
the results for fo^o, agree (within the large errors quoted there) 
with the ones in [27]. These large errors come mainly from 
mc,yu and (dd). The previous bounds can be used for excluding 
some experimental data and some theoretical estimates. 
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Table 7: Central values and corresponding errors for in units of MeV from 
the LSR and MSR at different values of the subtraction point in units of MeV 
for rhb = 7292 MeV. The +(resp. -) sign means that the values of fs increase 
(resp. decrease) when the input increases (resp. decreases). The Total error 
comes from a quadratic sum. 



II 

A* 


JBs 


t 


3 


mu (dd) (G^) 


Mq Ms (ss) Total 


LSR 
















225 


+ 18 


-1.9 +3.5 


-1.3 +1.4 +0.1 


-0.4 +0.4 +0.7 


19 


3 


226 


+ 17 


-8.5 +1.1 


-1.3 +1.1 +0.0 


-0.1 +0.5 +0.6 


19 


4 


230 


+ 10 


-8 +0.8 


-1.3 +1.1 +0.0 


0.0 +0.5 +2.9 


13 


5 


232 


+ 11 


-8.3 +0.7 


-1.9 +1.2 -0.3 


-0.5 +0.3 +2.5 


14 


6 


230 


+ 16 


-10.9 +0.8 


-0.9 +0.6 -0.2 


-OA +0.4 +2.5 


20 


7 


234 


+ 16 


-10.6 +0.7 


-1.2 +1.2 +0.1 


-0.3 +0.4 +2.7 


20 


MSR 














3 


224 


+8.0 


-14 +0.4 


-2.0 +1.4 -0.2 


-1.1 +1.0 +1.0 


16 


4 


235 


+ 13 


-10.6 +1.0 


-2.1 +1.3 -0.1 


-0.7 +1.0 +1.2 


17 


5 


240 


+ 12.4 


-19.2 +1.1 


-1.8 +1.2 0.0 


-0.3 +1.1 +0.7 


23 


6 


254 


+ 12 


-12.9 +1.0 


-2.0 +2.3 -0.8 


-2.4 +0.6 +1.5 


18 


7 


258.6 


+ 14 


-13.2 +1.0 


-2.4 +2.4 -0.5 


-2.2 +0.7 +1.5 


20 



Table 8: Results from the open charm and beauty systems in units of MeV and 
comparison with experimental data and lattice simulations using n/ = 2[63, 64] 
and rif = 3 [65, 66] dynamical quarks, fp are normalized as fj^ = 130.4 MeV 



Charm 


Bottom 


Ref. 




mbimb) 




1286(66) 


4236(69) 


This work 


1280(40) 


4290(140) 


ETMC[63] 


fo 


fB 




204(6) ^ 1.56(5)/, 


206(7) ^ 1.58(5)/, 


This work 


< 218.4(1.4) ^ 1.68(1)/, 


< 235.3(3.8) = 1.80(3)/, 


This work 


207(9) 




Data [2, 62] 


212(8) 


195(12) 


ET]V[C[63] 




193(10) 


ALPHA [64] 


207(4) 


190(13) 


HPQCD [65] 


219(11) 


197(9) 


FNAL [66] 




fs. 




246(6) ^ 1.59(5)A 


234(5) ^ 1.51(4)/, 


This work 


< 253.7(1.5) ^ 1.61(1)A 


< 251.3(5.5) ^ 1.61(4)A 


This work 


260(5.4) 




Data [2, 62] 


248(6) 


232(12) 


ETMC[63] 




219(12) 


ALPHA [64] 


248(2.5) 


225(4) 


HPQCD [65] 


260(11) 


242(10) 


FNAL [66] 



11. Summary and conclusions 

We have re-extracted the decay constants fD,D, and /b^b, and 
the running quark masses mc,b(^c,b) using QCD spectral sum 
rules (QSSR). We have used as inputs, the recent values of the 
QCD (non-)perturbative parameters given in Table 1 . The re- 
sults given in Eqs. (36), (39), (43) and (44) agree and improve 
existing QSSR results in the literature. Along the analysis, we 
have noticed that the values of the decay constants are very sen- 
sitive to the heavy quark mass and decrease when the heavy 
quark masses increase. Here we have used (for the first time) 
the scale independent Renormalization Group Invariant (RGI) 
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Figure 20: a) Upper bounds of fo^ (red triangle) and of fo (blue open circle) at different 
values of the subtraction point ji and for fnc-1467 MeV. The filled horizontal band is the 
average within the averaged error; b) the same as a) but for /^^ and fs with rhb -1292 MeV. 



heavy quark masses in the analysis. We have translated the on- 
shell mass expressions of the PT spectral function known to 
N2L0 into the MS one where (as has been already noticed in 
previous works [19]) the PT series converge faster. We have 
also remarked that fp and mq are aff'ected by the choice of the 
continuum threshold tc which gives the largest errors. Here, 
like in our previous works [13, 14, 19, 25], we have taken the 
conservative range of -values where the r- or ^-stability starts 
until the one where ones starts to have stability. We have also 
seen that the subtraction point ju aff'ects the truncated results 
within the OPE which has been the sources of apparent discrep- 
ancies and large errors of the results in the literature. Here, we 
have considered carefully the results at each subtraction point 
and deduced, from these "QSSR data", the final results which 
should be independent on this arbitrary choice. In view of pre- 
vious comments, we consider our results as improvements of 
the most recent ones to N2L0 and using MDA in [25-27]. 
The results on fo and fo^ agree within the errors with the data 
compiled in [2, 62], while the upper bound on fo^ can already 
exclude some existing data and theoretical estimates 
As one can see in Table 8, our results are comparable (in val- 
ues and precisions) with recent lattice simulations including dy- 
namical quarks [63-66] These agreements are not surprising 
as both methods start from the same observables (the two-point 
correlator though evaluated in two diflferent space-times) and 



^^Implications of the values of fD,Ds on the determinations of the CKM mix- 
ing angles have been discussed in details in [67]. 

^^A summary of the present results and comparisons with experiments and 
with the lattice ones will be presented elsewhere [1]. 



12 



use the 1st principles of QCD (here is the OPE in terms of 
the quarks and gluon condensates which semi-approximate the 
confinement regime). These agreements also confirm the ac- 
curacy of the MDA for describing the spectral function in the 
absence of a complete data, which has been tested earlier [5, 6] 
and in this paper from the charmonium and bottomium systems. 
MDA has been also successfully tested in the large A^^ limit of 
QCD in [36]. 
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